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Abstract 



Dynamics of spreading of small droplets on surfaces has been studied 
by the molecular dynamics method. Simulations have been performed for 
mixtures of solvent and dimer, and solvent and tetramer droplets. For sol- 
vent particles and dimers, layering occurs leading to stepped droplet shapes. 
For tetramers such shapes occur for relatively deep and strong surface po- 
tentials only. For wider and more shallow potentials, more rapid spreading 
and rounded droplet shapes occur. These results are in accordance with 
experimental data on small non - volatile polymer droplets. 

PACS numbers: 68.10Gw, 05.70.Ln, 61.20.Ja, 68.45Gd 



Experiments on the dynamics of spreading of tiny non - volatile liquid 
droplets have revealed a variety of droplet shapes on the molecular scale [|l], 0, 
|3|, £|] which cannot be explained by continuum hydrodynamic theories |5|, |, |7], 
ffl. These experiments have created a need for more microscopic models H. 
Some analytic models have been successful in describing dynamical layering 
|], [H], |TTJ. To obtain more microscopic understanding, molecular dynamics 
simulations have been performed fl^, [13[ . In particular, Nieminen et al. [|l4f 
simulated Lennard - Jones (LJ) particles and observed the appearance of a 
precursor film, dynamical layering, and a crossover from almost linear (~ t a9 ) 
to t ' 5 behavior for the width of the precursor |7|, ||. 

However, recent experiments using a variety of polymer liquids on en- 
ergetically different surfaces show features which are beyond the simple LJ 
model of a liquid |], 0, |3|, [|. For example, a droplet of squalane develops only 
a tiny foot preceding the droplet 0] , whereas in the case of polydimethylsilox- 
ane (PDMS) a monolayer thick precursor layer is seen [0], except on silver 
where the droplet has a spherical cap and Gaussian tails at late times H]. 
Furthermore, tetrakis (2 - ethylhexoxy) - silane develops multiple (up to four) 
precursor layers and a tiny cap in the center of the droplet |TJ . Venalainen et 
al. Jl5j have recently demonstrated that most droplet shapes can be qualita- 
tively reproduced by an effective solid on solid model where interparticle and 
substrate interactions are varied. However, there still remains an obvious 
need to understand the role of microscopic solid - liquid and liquid - liquid 
interactions, temperature, and, in particular, the chainlike morphology of 
the molecules. Although most experiments on tiny droplets have been done 



with rather short molecules beneath the static entanglement regime [16], re- 
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cent work shows || how spherical droplet shapes can occur due to strong 
blocking effects between short chainlike molecules. This is the motivation for 
the present work, where we present a study where the chainlike structure of 
the spreading liquid and the influence of the surface potential are taken into 
account microscopically within the molecular dynamics (MD) method. 

We consider a model system of a mixture of "solvent" (single atom) parti- 
cles and "polymers" (short effective chainlike molecules) which interact with 
each other via an LJ potential, with potential parameters 07 and e/ for the 
width and depth of the potential, respectively. The substrate is modeled by 
a flat continuum LJ material ||14|| . For the substrate - fluid interaction, the 
depth of the potential is comparable to e s a^ and the width to a s (see Ref. 
fLj). The n - mers in the system consist of n LJ particles, which are in- 
terconnected by a very rigid but orientationally isotropic harmonic oscillator 
pair potential V c = \k(r — r ) 2 , where k = lOOe f/crj. There is also an angle 
dependent potential Vg = eg (cos 9 + 1), which favors linear molecules. Here 
we present results for eg = 10ej which turns out to be the most interesting 
case; results for more flexible molecules are similar, and will be published 
elsewhere PT7| . Due to the strong pair potential between subsequent atoms 
comprising an n - mer, a strong repulsion is created if any n - mers try to 
spatially overlap. 

The dynamics of the system is described by the usual equations of mo- 
tion |T4}1 : dri/dt = Pi/rrii, dpi/dt = ViV — wpi, and dw/dt = [Sj(p^/mj) — 
NkT]/NkTr 2 , where iV is the number of degrees of freedom, T the tempera- 
ture for the thermostat, and r is a relaxation time. The time scales are here 
chosen in the same way as in Ref. Initially a ridge - shaped droplet 
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is constructed with periodic boundary conditions along the direction of the 
ridge. The spreading takes place in the direction perpendicular to the ridge. 
The equations of motion are solved using modified velocity Verlet algorithm 
(see e.g. @ and JOJ). 

The initial configuration is constructed so that the n - mers are bent 
90 degrees at each joint, while their directions are arbitrary. To attain an 
optimum packing density the system is compressed, after which the system 
is allowed to find a lower energy state. After finding the best possible con- 
figuration we switch on the substrate interaction and allow the system to 
evolve. 



In the previous study of Nieminen et al. [|14| the temperature was kT s 



0.8e/, which is well above the triple point of an LJ material pOfl . In the 
present case the larger mass and binding energy of the n - mers allows us to 
systematically vary both the surface potential and the temperature, and thus 
e.g. study the effects of changing the viscosity of the liquid. In the present 
work, we employ several sets of substrate interaction parameters: e s = 5e/ 
and a s = <jf (used in Ref. fl^l), denoted as Sr, e s = e/ and cr s = 5.O07 
(Sn)'i e s = 0.02e/ and a s = 7.3o/ (Sjjj); and e s = 0.002e/ and a s = 7.3<T/ 
(Siv)- The differences between these potentials can be characterized by their 
steepness (gradient) (~ e s a s 2 ), depth, and width. The potential Sji is much 
steeper and deeper than Sj, while Sjh and SVv are both relatively shallow, 
the latter being particularly flat. 

First we describe results for mixtures of solvent and dimer particles. Qual- 
itatively, a dimer concentration of 50% or 100% does not alter the essential 
features of spreading as compared to a pure solvent droplet. Dimer droplets 



3 



exhibit well separated layers such as observed for solvent droplets for Si, 
Sjj and Sjjj. However, in order to clearly observe this, a somewhat deeper 
substrate potential Su must be used (for Sj only one precursor layer forms). 
Despite this slight difference dimer droplets show a rather clear crossover 
from the almost linear to t 5 spreading as described in Ref. |TJ| for solvent 
droplets. The presence of dimers, however, makes the spreading dynamics 
slightly slower. 

In the case of tetramers, effects arising from the chain structure of the 
liquid become more apparent. The first striking feature is a dynamical or- 
dering of the rather stiff chains along the "flow" of the spreading droplet, 
perpendicular to the surface (Fig. 1). The nature of layer separation is also 
changed. Only for the very strong surface potential Su a precursor layer 
with thickness of about one monolayer is formed. When surface potential 
becomes more shallow, a distinct change in the droplet shapes occur. For 
Sin a vertically continuous precursor layer whose effective thickness is two 
layers is seen. When the substrate potential is further weakened to Sjv, the 
two layers tend to merge together. This leads to a noticeably more rounded 
droplet configurations than for dimers. The rate of spreading of the droplet 
also becomes faster from Su to Sui- This is consistent with the observation 
that for a strongly attractive, narrow surface potential the tetramers become 
tightly packed parallel to the surface, which slows down their migration in 
the effectively two dimensional layer ||. Such layers also remain relatively 
well separated. On the other hand, a wider surface potential does not force 
the tetramers in layers, thus allowing faster dynamics and formation of a 
more rounded droplet configuration. Mixing of the layers in Fig. 1 is seen 
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which leads to a very complex structure within the precursor layer. When 
the substrate interaction is still weakened from (Sin) to (S/y), no further 
increase in the rate of spreading is not seen. 

We have also studied the effect of increasing tetramer concentration start- 
ing from a pure solvent droplet, as shown in Fig. 2. This slows down the 
spreading rate in analogy with the dimer case. Most importantly, we find that 
the width of the precursor film w obeys a finite size scaling form as a function 
of the number of LJ particles N 0, |17| w = t x (j)(t/N y ), where x « 9/10 0, 
and y ~ 0.67 jx. This was checked for 2380 < N < 4080. This indicates that 
the spreading is not purely diffusive but rather caused by a combination of 



the pressure of the cap and collective migration of the tetramers ||21|| . The 
scaling function <p(z) ~ const, for z < 1, and 4>(z) ~ z 1 / 2 -* f or z ^> \ 
Details of these results will be published elsewhere Jl7| . 



Next we study the density profiles of the droplets in order to make com- 
parison with experiments 0, |3], |]. As seen in Fig. 3(a), the density profile 
for a solvent droplet shows a stepped precursor layer for Sn, indicating sepa- 
ration of layers. Results for dimers are analogous, with slightly more rounded 
step edges for weaker potentials. For tetramers in the case of Sn, the density 
profile shows a stepped shape, indicating dynamical layering. For weaker po- 
tentials the comparison of snapshots and profiles shows mixing of the layers. 
For Shi th e precursor film on the whole is fairly thick and the density con- 
tinuosly increases towards the centre. The density profile is seen to develop 
towards a rounded cap JJ, which is most clearly seen in Fig. 3(b) for the case 
of Siv- Additional studies with more flexible chains and octamers support 



these conclusions [17 
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The profiles of Fig. 3 bear a remarkable qualitative resemblence to some 
of the experimental ones. The stepped profile of Fig. 3(a) is qualitatively 
similar to that of a tetrakis droplet consisting of star - shaped molecules 
Our results indeed suggest that for spherical or dimer molecules separation 
of layers tends to occur. On the other hand, the profile of the tetramer 
droplet (Fig. 3(b)) resembles the profiles of a PDMS droplet on silver ||, 
and on silicon ||. These more rounded shapes should occur for chainlike 
molecules, when the surface attraction is not strong. Furthermore, on sil- 
icon a change of the droplet from a rounded towards a stepped shape was 
observed on a "high - energy" surface. Our results suggest that this indi- 
cates increased substrate attraction, which tends to separate the layers. A 
quantitative comparison of our results with experiments is difficult, however, 
since experimental droplets are microscopic in the vertical direction only. 
Nevertheless, our results support the conclusion that weakening the surface 
potential causes faster spreading and more rounded droplet shapes for poly- 
meric liquids. Faster spreading of the precursor film for weaker potentials 
has also been predicted in Ref. 0. 

Finally, we have also studied the effect of temperature. As expected, 
increasing the temperature increases the spreading rate. For solvent droplets 
{Sjj) changing the temperature from T = T$ to T = 1.5Ts simply makes the 
steps more rounded. For tetramer droplets in the case of Sni, the effects are 
relatively small. For T = 4Ts the profile shows a minimum at the centre of 
the droplet, after the vanishing of the cap. This is followed by a fairly flat 
rounded profile, as seen in Ref. ||. 

To summarize, in this work we have studied effects arising from the chain- 



6 



like molecular structure within spreading droplets, and the form of the surface 
potential. Our model simulations suggest that for relatively strongly attract- 
ing surfaces, dynamical layering and stepped droplet shapes tend to occur. 
On the other hand, the mixing of chainlike molecules opposes this trend lead- 
ing to more rounded droplet shapes on weaker substrates. Our results also 
support the conclusion that weaker substrate attraction enhances the rate of 
spreading. These results are in agreement with experiments Jl], j|, [3|, |3J]. We 
hope that our work inspires further systematic experiments on the effects of 
chain lengths and surface potentials on the spreading of tiny droplets. 

Acknowledgements: We wish to thank K. Kaski, R. Swendsen, and O. 
Venalainen for useful discussions, and S. Herminghaus for a critical reading 
of the manuscript. This work has been supported by the Academy of Finland. 
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Figure Captions 



Fig. l(a)-(b). Snapshots of a projection of a three dimensional pure tetramer 
droplet for potential Sjv, with N = 2380. Times correspond to 50 and 75 in 
reduced units ||14j| . Mixing of the layers is clearly visible. Units in this and 



the other figures have been obtained from aj = 2.6 A. Due to the effective 
nature of the chainlike molecules the length scales should be considered ef- 
fective only. 



Fig. 2. The effect of tetramer concentration c on the time evolution of the 
precursor width w for potential Su' c = (full line), c = 0.5 (dash-dotted 
line), and c = 1 (dotted line). 

Fig. 3. Smoothed density profiles for two different droplets, (a) A solvent 
droplet for Su and (b) a tetramer droplet for Siy. The relative times are ti, 
4ti, and 6ti. 
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